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Abstract 

When a solid body is freely rotating at an angular velocity fi, the ellipsoid 
of constant angular momentum, in the space f^i, Vl2-i ^3; has poles corre- 
sponding to spinning about the minimal-inertia and maximal- inertia axes. 
The first pole may be considered stable if we neglect the inner dissipation, 
but becomes unstable if the dissipation is taken into account. This happens 
because the bodies dissipate energy when they rotate about any axis differ- 
ent from principal. In the case of an oblate symmetrical body, the angular 
velocity describes a circular cone about the vector of (conserved) angular 
momentum. In the course of relaxation, the angle of this cone decreases, 
so that both the angular velocity and the maximal-inertia axis of the body 
align along the angular momentum. The generic case of an asymmetric 
body is far more involved. Even the symmetrical prolate body exhibits a 
sophisticated behaviour, because an infinitesimally small deviation of the 
body's shape from a rotational symmetry (i.e., a small difference between 
the largest and second largest moments of inertia) yields libration: the pre- 
cession trajectory is not a circle but an ellipse. In this article we show 
that often the most effective internal dissipation takes place at twice the 
frequency of the body's precession. Applications to precessing asteroids, 
cosmic-dust alignment, and rotating satellites are discussed. 
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I. INTRODUCTION 



A complex rotational motion of a free solid body is an evidence of its rotation-axis' 
wobbling about the angular momentum. Indeed, a solid body in a long-established regime 
of free rotation must have its axis of rotation parallel to the angular momentum: this 
configuration will minimise the kinetic energy, the angular momentum being fixed. The 
body achieves this minimisation by aligning both its axis of rotation and axis of maximal 
inertia parallel to the angular momentum. By the end of this relaxation the body comes to 
steady spinning about its maximal-inertia axis. Any deviation from this regime witnesses 
either of the influence of the tidal forces, or (in the case of comets) of the result of jetting, 
or of an impact experienced by the body within its characteristic time of relaxation, or of 
the entire body being a wobbling fragment of an asteroid disrupted by a collision (Giblin 
& Farinella 1997 )\ (Giblin et al 1998)2, (Asphaug & Scheeres 1999)^. The contest 
between the impacts (or the tidal forces, or the cometary jetting) on the one hand, and 
the relaxation mechanism(s) on the other hand, determines the dynamics of the body 
rotation. 

A study of the rotation of asteroids and comets may thus provide valuable infor- 
mation about their recent history. Several examples of complex motion have already 
been registered. Among the asteroids, 4179 Toutatis furnishes another example of 
wobble (Ostro et al 1993)^ (Ostro et al 1995)^ (Ostro et al 1999)^ see also 
|http: / / www.eecs.wsu.edu^hudson/ asteroids. html| ). Among the comets, P/Halley is cer- 
tainly an example of such a tumbling object (Sagdeev et al 1989)^, (Peale & Lissauer 
1989)^ (Peale 1991)^ (Behon et al 1991 (Samarasinha & A'Hearn 1995)". An- 
other example is comet 46/P Wirtanen. ROSETTA mission is supposed to explore this 
comet soon (Hubert & Schwehm 1991)^^. The spacecraft will be carrying a high-resolution 
OSIRIS imaging system (Thomas et al 1998)^^ that, probably, will enable the astronomers 
not only to see the precession but also to observe its relaxation during the 2.5 years of 
the spacecraft's escorting the comet. Our estimates^^'^^ show that the angular resolution 
of the currently available equipment gives to such sort of experiment a good chance of 
success, provided the experiment includes measurements performed at least half a year 
apart from one another. An important factor influencing the success or failure of such an 
experiment is the jetting intensity of the particular comet. 

Another field of application for this study is the cosmic-dust alignment: some of the 
alignment mechanisms are very sensitive to coupling between the angular velocity and 
angular momentum of the interstellar grains (Lazarian & Draine 1999)^^, (Lazarian & 
Efroimsky 1999) 

The third possible application of the developed formalism could be spin stabilisation of 
spacecrafts, including spacecrafts with a precession damper (Chinnery and Hall 1995)^^, 
(Hughes 1986)^^, (Levi 1989)^^. An interest in studies of nonrigid-body dynamics with 
applications to spacecraft motion emerged after launch of "Explorer" satellite in 1958. 
(I am thankful to Vladislav Sidorenko for drawing my attention to this example.) The 
satellite was a very prolate body with 4 small deformable antennas on it. It had been 
supposed that it would rotate about its minimal-inertia axis. Instability of this motion 
was a major surprise for mission experts (Modi 1974)^''. 

On general grounds, the necessity of relaxation is evident: the system must reduce its 



2 



FIGURES 

03 




FIG. 1. The ellipsoid of eonstant angular momentum in the angular- velocity space 
(f^i , $72) ^s)- Lines on its surface denote its intersections with the ellipsoids appropriate 
to different values of the rotation energy T^in- Pole A corresponds to the maximal value of 
Tkin and, thereby, to rotation about the minimal-inertia axis. Pole C corresponds to a complete 
relaxation, when T^in is minimal and the body is spinning about its major-inertia axis. In the 
course of precession, the angular-velocity vector ft moves along the constant-energy lines, and 
slowly shifts from one line to another, in the direction from A to C. The above picture describes 
the case of an almost prolate body: ^ I2 > Ii- The trajectories are circular near A and, 
in this case, remain almost circular all the way from A to the separatrix. Crossing over the 
separatrix may yield chaotic fiipovers, whereafter the body will begin libration. The trajectories 
will again become circular only in the closemost vicinity of C. 

kinetic energy down to the value that is minimal available for a fixed angular momentum. 
What particular physical effects provide this relaxation? One phenomenon, relevant to 
tiny grains (like those of the interstellar dust) but feeble for large samples, is the Barnett 
dissipation called into being by the oscillating (due to the precession) remagnetisation of 
the material, caused by the Barnett effect. (See, for example, (Lazarian 1994)^^, (Lazarian 
& Draine 1997)^^, and references therein.) Another process, relevant in small granules, 
and overwhelmingly leading for large bodies, is the inelastic dissipation. It is produced 
by the precession-caused alternating stresses and strains. 

A pioneering paper on the inelastic dissipation was published by Burns & Safronov^^ 
back in 1973. Later the inelastic dissipation in small and in large freely-rotating oblate 
bodies was addressed in (Lazarian & Efroimsky 1999)^^ and (Efroimsky & Lazarian 
2000)^^, correspondingly. In the present article wc shall tackle to the dynamics of a 
body of arbitrary values of its moments of inertia. The issue is nontrivial. For example, 
a perfectly biaxial prolate spheroid has two axes of equal maximum moment of inertia 
and, therefore, really does not have a stable rotation pole. For a triaxial figure of a shape 
slightly deviating from a symmetrical prolateness, the situation is that the precession tra- 
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jectory is not a circle but an elongated ellipse, with long axis around the "waist" of the 
"cigar" figure (Fig.l). For the tumbhng case, the spin axis circulates all the way around 
the body, and in fact the circulation is more or less around the long axis. This leads to 
the curious effect, observed in computer simulations of the rotation of asteroid (433) Eros, 
that the body appears much of the time to be spinning nearly about its long axis (Black, 
Nicholson, Bottke, Burns & Harris 1999)^^. 

In what follows we briefiy review the main facts and formulae describing the solid-body 
rotation (section II), dwelling comprehensively upon the case of an almost symmetrical 
prolate body (section III). We divide the motion into four distinct stages (section IV). 
Then we discuss the relaxation rate (section V) and explain the nature of the nonlin- 
earity emerging in this problem (section VI), whereafter we compute the stresses arising 
in a processing body, and calculate the energy density of the appropriate deformations 
(section VII). In sections VIII and IX we calculate the rate of internal dissipation. In 
section X wc draw conclusions and mention some practical applications of the formalism 
developed. In section XI we briefly account of the vississitudes of the generic case. 

II. NOTATIONS AND ASSUMPTIONS 

We shall discuss free rotation of a solid body, using two Cartesian coordinate systems, 
each with an origin at the centre of mass of the body. The inertial coordinate system 
(X, Y, Z), with unit vectors ex, ey, ez, will have its Z axis parallel to the (conserved) 
angular momentum J. Coordinates with respect to this frame are denoted by the same 
capital letters: X, Y, and Z. We shall also use the body frame defined by the three 
principal axes of inertia: 1, 2, and 3, with coordinates x, y, z and unit vectors ei, 62, 63. 

We denote the angular velocity by f2, while u; will stand for the rate of precession. 
Components of f2 in the body frame will be called ili,2.3 • 

Due to the lack of an established convention on notations, we present a table that 
hopefully will prevent misunderstandings. 
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TABLES 





Principal 
moments 
of inertia 


Components of 

angular 

velocity 


Frequency of 

precession 


(Purcell 1979), 
(Lazarian & 
Efroimsky 1999), 
(Efroimsky 2000), 
(Efroimsky & 
Lazarian 2000), 
present article 


^3 > > h 






(Synge & Grif- 
fiths 1959) 


A>B>C 


wi , a;2 , 


P 


(Black et al. 
1999) 


C>B>A 




V 
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A free rotation of a body obeys the Euler equations 

^ (Ji a) = (/,■ - 4) Qk (2.1) 

where {i,j,k) = (1,2,3) (with cychc transpositions), and the principal moments of 
inertia are assumed to obey: 

h < h < h . (2.2) 
In the approximation of an absolutely solid body, the equations simplify: 

li hi = {ij - h) fife . (2.3) 

This neglect of li Qi against Jj Qi does need a justification, because the inelastic 
relaxation we are going to describe is due to small deformations that yield nonzero J, . 
To validate the neglection, i.e., to prove that li/Ii <^ Cli/Qi , one must recall that for a 
rotational period r 

- r-' , ^ r-'e (2.4) 

e being a typical value of the relative strain (which in real life rarely exceeds 10~^). 
Conservation of the angular momentum J and the kinetic energy Tfcj„ entails: 

I^Ql + Ilnl + Il^l = J2 , (2.5) 



h nl + hnl + hnl = 2 (2.6) 

In the context of our study (which is aimed at estimating the rate of relaxation), equation 



(2.6) is applicable as long as we accept the adiabatic approach, i.e., assume the relaxation 
to be a "slow process", compared to the "fast processes" of rotation and precession. A 
rigorous formulation of this assertion is based on formulae ( |8.24| ) and ( |9.15|) obtained 



below. These formulae are to become the main result of our article. They give the 
relaxation rate d< sin^ 6 >/dt , 9 being the angle between the angular momentum and 
the major-inertia axis of the body, and averaging being performed over the precession 
period (Appendix A, formula (Al)). The exact formulation of the adiabatic approach 
will read: 

d < sin^ 6 > , 

It — « ^ • (") 

uj being the precession rate. This condition will adumbrate the applicability realm of the 
solutions ( p.24| ) and ( |9.15| ) to be derived. 

Equations ( p.5|) and (plGp may be resolved with respect to and ^2 • 



H V\ — H) h ^1 — ^3 

n2 D cn2 D — '^hTkin - Q _ h h - h 

- K - ti\l2 , n = — - , ^ = — 

-'1 l-'s ~ -'ij -'1 -'3 ~ -'1 



(2.^ 
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substitution whereof in ( p.3|) , for ^ = 2 , gives: 



hi 



p - Qnfj (i? - snl) 



(2.9) 



It is possible (Synge & Griffith 1959)^^ to pick up such positive functions P , u , k 
of the arguments /i,2.3, Tkm, that the rescaled time and second component of the 
angular velocity, 



t' = OOt , ^ = ^2/(3 



satisfy 



(d^/dtf 



Oil - ee 



(2.10) 



(2.11) 



with k < 1 . Solution to this equation is the Jacobian elliptic function sn(t') , so that 



(3 sn uj{t — to)? 



(2.12) 



to being an arbitrary constant. It is known (Ibid.) that substitution of the latter in 
yields, for 3"^ > 2 h Tk 



■ kin 



7 cn 



uj{t — to), k^ 



a dn 



uj{t — to)? k^ 



(2.13) 



(Recall that what we call VL^ , in (Synge & Griffiths)^^ is called cui, while (Black et al) 



25 



denote it Uc-) In the above formulae 



a 



J2 - 2 Ji T; 



kin 



\ h{h-h] 



ijj 



(J2 - 2 h Tuin) ih - h) 



V h h h 

while for < 2 I2 T^in one arrives to: 

ill = J dn i^{t — to), k^ 

where 



2 ^3 Tkin — J" 
\ hih-h) 

k = ^ 
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2 -^3 Tfcjn 
\ hih-h) 



II2 - h 2/3 T; 



kin 



I2 — 2 Ii Tkin 



a cn 



u{t — to), k^ 



(2.14) 



(2.15) 



a 



2/iT; 



kin 



1(2 hh 



P 



\ hih 
ih - h) 



h h h 



2hJ) 

h] 

k 



kin 



7 



2 h T, 



3 kin 



\ hih-h] 



I2 P 

h 2hh 



2hh 



kin 



(2.16) 



kin 



In some books (like for example in Abramovitz & Stegun 1964)^'' notation m = k^ is 
used 
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Mind that in Synge & Griffith'^^ expression (14.116a) for 7 is given with a wrong 
sign. Our expression for 7, as given by our formula (|2.14| ), makes the expressions (|2.12 - 
|2.13| ) coincide, in the hmit of oblate symmetry ( Ji = I2), with the well-known Eulerian 
solution: Qi = Q± cosuj{t — to) , ^2 = ^± sinuj{t — to) , ^3 = const , for > . 
Our choice of signs is correct since it leaves f2 parallel to J in the relaxation limit. 

Our ultimate goal is to compute the rate of inelastic dissipation caused by alternating 
stresses in a wobbling body. To know the picture of stresses, one should begin with 
derivation of the acceleration experienced by a point (x, y, z) inside the body. We mean, 
of course, the acceleration with respect to the inertial coordinate system (X, Y, Z) , 
but for convenience of the further calculations we shall express it in terms of coordinates 
{x, y, z) of the body frame (1, 2, 3) . The position, velocity and acceleration in 
the inertial frame will be denoted as r, v, a , while those relative to the body frame 
(1, 2, 3) will be called r', v' and a' (where r' = r ). The proper acceleration (i.e., 
that relative to the inertial frame) will read: 



a' + ri X r' + 2 ri X v' 



n X (fix r') 



(2.17) 



In the beginning of this section we justified, on the grounds of the strains being small, 
our neglect of Ij Q against Ij Cl . In a. similar manner we shall justify the neglection 
of the first and third terms in the right-hand side of the above formula: rotation with a 
period r , of a body of size / will yield deformations 61 ^ el and deformation- caused 



61 /t ^ el/r and accelerations 



61 /t^ = el/r^ , e being the 



velocities v' 

relative strain. We see that v' and a' are much less than the velocities and accelerations 
of the body as a whole (that are about 1/t and //r^ , correspondingly). Henceforth, 

a ^ 17 X r' + f2 X (f2 X r') = 

ei { (I2 z - (l^y + ^2 i^iV ~ ^2 x) - f23 (^3 x - ^li z) } 
+ 62 { 0,3 X — Qi z + Q3 {Q2 z — Q3 y) — Qi {Qi y — Q2 x) } 
+ { Qi y — Q2 X + Qi {Q3 X — Qi z) — Q2 (^2 z — Q3 y) } = 

= ei { tl2 z - (is y - X (ill + ilj^ + y ill il2 + z ill n^} 

+ e2 {tl3 X - ill z - y (ill + ilj^ + Z ils il2 + X ill il2} 

+ e3 {iliy - il2X - z (ill + ill) + xiliil3 + yil2il3} ■ (2.18) 



where, according to ( p.8| ) and ( |2.12| ): 

ill + ill = P+i'^-Q) f^^ sn^ Wit - to), P 



(2.19) 



ili + ill = R + (1 - S) sn' iuit-to), k 



(2.20) 



+ (]f = (P + i?) - (g + S) (3' sn' uj{t - to), k 



(2.21) 
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^3^1 = a 'J dn uj{t — to)) cn uj{t — to), k 



CI2 = P uj cn uj{t — to), dn Lu{t — to), k 



the derivative being written in compliance with (Abramovitz & Stegun 1964) 
16.16.1. The other items emerging in (|1|) will read, for > 2 Is 



^3 = O! P dn uj{t — to), k"^ sn uj{t — to), k 



0.2^1 = P 'J cn uj{t — to), k'^ sn uj{t — to), k 



(li = — 'J CO sn uj{t — to), k'^ dn uj{t — to), k 



and 



VL-x = — a uj k cn 



uj{t — to), k'^ sn uj{t — to), k 



while for < 2 12 T^m ■ 

fls Q2 = C( P cn uj{t — to), k"^ 



sn 



uj{t — to), k^ 



VL2VL1 = P 'J dn uj{t — to), k"^ sn uj{t — to), k 



Vti = — J uj k'^ sn uj(t — to), k'^ cn u{t — to), k 



and 



a LO sn 



^{t — to), k'^ dn uy{t — to), k 
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III. ALMOST PROLATE BODY 



Before pursuing the generic case, let us dwell for a minute on the case of an almost 
symmetric prolate top of and I2 having close values: 



(3.1) 



I.e., 



h ~ h 



h » /3 



Jo ~ 



(3.2) 



As we saw above, the solution depends upon the sign of ( — 2 12 Ttin) ■ According to 
(O-O), 



kin 



^3 -^3 ~ ^3 -^3 -^2 + if — h I2 — ^3 li 



J2 _ ^hh 



(3.3) 



Based on (3.1 - 3.2), one can introduce the following parameters: 



(3.4) 



and 



h 



(3.5) 



Besides, we shall use a time-dependent quantity ^f/^l ■ It is not, of course, a geometrical 
parameter worthy of the name, though formally one may consider it as a sort of parameter 
in ( p.3| ). This wanna-be parameter may be used as a measure of the system's approaching 
the steady regime: if Ql/fl'^ -C 1 over the entire time of one wobble then the relaxation 
is almost over, and the angular velocity fl is precessing about J , with a small amplitude 
(i.e., describing a narrow cone). Below, the meaning of the words like small and narrow 
wiU become understandable. Meanwhile, one can write down (13. 31) as 



2 I2 Tkin 



fig J3 



fi2 



(3.6) 



and easily find that for fixed values of q and s (i.e., for a particular prolate body) the 
narrowness of the precession cone yields: 



2/2T; 



kin 



nl il q 



This approximation becomes true at the late stage of relaxation, when 



1 



.(1 



(3.7) 



(3. 
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holds through the duration of one wobble. Assume, following (Black, Nicholson, Bottke, 
Burns & Harris 1999)^^, that the moments of inertia of asteroid (433) Eros relate as 
Ii : I2 : I3 = 1:3: 3.05 . With these numbers plugged in, the above formula 
will give: l^i/fisl <ti 0.4 . Formally, ( p.7| ) was derived from ( |3.6| ) by keeping q and s 
fixed, and making the "parameter" Ql/Qf approach zero. After this is done, one may 
consider a variety of geometries, and make q approach zero. Then — 2 J2 Tj^in will 
approach zero, always remaining positive, so that the end of relaxation will be described 
by the solution (|2.12| , |2.13| ) with frequency u expressed by ( |2.14| ). 



On the other hand, one might as well perform a different, unphysical trick: for a 
fixed fli , begin with a limit g ^ , and only afterwards choose the case of the 
small-amplitude wobble (i.e., consider fif/fig — > ). The first limit will give: 



r - 2h ^ -nii'i^sii-s) , (3.9) 

"3 



After the second limit is taken, 3^ — 2 I2 Tkm will approach zero, always remaining 
negative, so that the final stage of relaxation would obey the solution (|2.12| , |2.15|) with 
frequency uj expressed by ( ^.161 ). 



We see that the operations g — > and ^ do not commute. From the 

physical point of view, an observer studying, for a variety of samples, the end of relaxation 
should first fix the shape of the body (i.e., assume, for example, that q is small enough 
but constant). Only afterwards he may state that he is interested only in the final spin, 
i.e., assume -C q s^^ (1 — s — q)^^ . As explained above, this observer will see 

that the end of relaxation takes place at frequency uj expressed by ( p.l4|) . 



An opposite order of limits would be physically meaningless, in that it would not help 
us to describe the behaviour of a particular body. 

We had to dwell on this issue so comprehensively because it would be very important 
to understand better the following two statements made in (Black, Nicholson, Bottke, 
Burns & Harris 1999)2^: 

1) . In the limit where Ii = I2 , the circulation region around the maximal-inertia 
axis vanishes, and 

2) . All trajectories circulate around the minimal-inertia axis. (This statement is 
fortified by the following argument: "the shghtest perturbation would cause such an 
object to "roll" about its long axis.") 

To analyse this statement by Black, Nicholson, Bottke, Burns & Harris let us cast 
it in a more exact form. First of all, we should understand which of the two possible 
sequencies of limits these authors implied. In case their statement implied the limit 
g — >■ taken first, and the relaxation limit Qf/Ql — > taken afterwards, then 
according to ( p.l6|) the frequency of precession u will read as 



'(2 Is Tkin - J2) ih - h) K^l h ih - h) + h ih - h)) ih - h] 



h h h V h h h 



\ — T~~7~~ ~ l^il ' f'^^ g — > and being fixed. 

V h h 

It will then approach zero as . However, as explained above, such a sequence of 

limits would be unphysical, while a physical way is to fix the body shape first (i.e., to fix 
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the difference {I3 — I2) ), then to take the relaxation hmit —>■ , and only after 



that to let q approach zero. In this case, ( |2.14| ) will yield u approach zero as y^g , i.e., 
as - h ■ 



h ih - h) + ^^i h ih - h)) ih - h 



h h h 



I a. 



for Vt2 



h h h 



and the shape being fixed. 



The above, physically meaningful, expression coincides with formula (2) in (Black, Nichol- 
son, Bottke, Burns & Harris 1999)^^, which means that the authors chose the right se- 
quence of limits. We certainly agree with the first of the above two statements derived 
by the authors from this formula: in the limit of 1^ = 12, the circulation region around 
the maximal- inertia axis vanishes. To understand the second of the above two statements 
made in (Black, Nicholson, Bottke, Burns & Harris 1999)^^, i.e., that "the slightest pertur- 
bation would cause such an object to "roll" about its long axis" , note that the relaxation 
limit was achieved in (|3.11| ) by letting ^2 vanish, with no assumptions made about Qi . 
This happened because u remains fii— independent as long as ( |2.14| ) (and therefore 
( |3.11|) ) may be used. These formulae may be used when the right-hand side of (|3.3|) (and 
that of (P^)) is positive, i.e., when in ( |3.8| ) we have at least "<" (not necessarily "<^"). 
In the case of asteroid (433) Eros, for example, this will work as long as, approximately, 
< 0.4|f23|. Thus we must agree with the statement about "rolling" around the 
minimal-inertia axis, but we have to add an important comment to it: 

As long as this rolling is slow enough, it will leave the spinning body within the realm 
of solution ( 2.13^ ), ( 2.14 ) appropriate to the final stage of relaxation. Too fast rolling will 
make it obey a different solution, ( ^.ISj ), with oj expressed by ( ^.Idj ). 

As an illustration, let us consider, in the space Qi , Q2 , ^3 , the angular-momentum 
ellipsoid 



I2 ^2 



+ 



Jg fig 



(3.12) 



and mark on its surface the lines of its intersection with the kinetic-energy ellipsoids 



2 Tkin — Ii fll + I2 ^2 



Jg fig , 



(3.13) 



at different values of energy. Let, on Fig. 1 , the starting point of motion be somewhere 
close to the pole A : the vector f2 initially is almost perpendicular to the major-inertia 
axis (3). Precession of the body will correspond to vector f2 describing a constant- 
energy "circle" in Fig. 1 . The word "circle" is standing in quotation marks because 
this trajectory is circular as long as we do not approach the separatrix too close. In the 
cause of the body rotation, f2 will be gradually changing the "circles" it describes, and 
will eventually approach the separatrix, en route whereto the "circles" will be getting 
more and more distorted. If the dissipation is slow, i.e., if the kinetic-energy loss through 
one precession period is less than a typical energy of an occational interaction (like, say, 
a tidal-force-caused perturbation) then chaotic motion is possible when f2 is crossing 
the separatrix: the body may perform flip-overs. After that the body will embark on the 
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stage of tumbling. As one can see from Fig. 1 , the tumbling will eventually turn into 
the final spin, i.e., into an almost circular small-amplitude precession of f2 around point 
C . However, this point will never be reached because the alignment of f2 along J has 
a vanishing rate for small residual angles: it is evident from formulae ( p.24| ) and ( |9.15| ) 
below, that at the end of the relaxation process the relaxation rate approaches zero, so 
that small-angle nutations can persist for long times. 

Now we understand that if the afore mentioned "rolling" becomes too swift, this will 
look as a "jump" over the separatrix in Fig.l. If we assume that {I3 — /2)//3 is 
infinitesimally small, then the separatrix will approach pole C infinitesimally close, and 
the smallest tidal interaction will be able to push the vector Q across the separatrix. 
In other words, the body, during the most part of its history, will be precessing 
about its minimal-inertia axis. For the first time this fact was pointed out in (Black, 
Nicholson, Bottke, Burns & Harris 1999)^^. 

Dependent upon the particular value of (Is — l2)/h and upon the intensity of the 
occational tidal interaction, the vector Q will be either driven from pole C back to 
the separatrix, without crossing it, or will be forced to "jump" over it. In the latter case, 
chaotic flipovers will emerge while fl is crossing the separatrix. 

As already mentioned, f2 will never approach pole C too close because in the 
vicinity of C the relaxation rate asymptotically vanishes. The behaviour of fl after 
its crossing the separatrix will be determined by two factors: on the one hand, occasional 
tidal interactions will push Q, towards or over the separatrix; on the other hand, the 
inelastic-dissipation process will always drive fl in the direction of pole C (though, 
once again, it will never manage to bring it too close to C , for the above mentioned 
reason). Sometimes this regime will be interrupted by collisions which may drive fl far 
away from the separatrix, in the direction of pole A . 



IV. STAGES OF MOTION 

We shall consider motion of a freely rotating body moving through four stages of 
relaxation. 

The first stage will be called "the initial spin" . It begins when the body rotates about 
some axis (almost) perpendicular to that of major inertia. This motion is characterised 
by negative (J^ — 2 I2 Tkin) , so that the frequency oj and parameter k are expressed 
by (|2.16|) . The initial spin starts when Vl\ is small. Namely, according to (|3.3|) , 



^ > ^3-/2 h ^ q_ 

h-h h s{l- q 



2 > = \ .^ • (4-1) 



Note that does not enter this condition at all. 

Now suppose that the body starts its rotation about an axis that is, for example, close 
to the minimal-inertia axis, so that is much less than Vl\ . Then, according to (|2.5|) , 
(PD and (|2l^ ). 
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(1 - s) 



(4.2) 



and 



h - h h ih - h) + h ih - h) 

the inequality ensuing from ( |3.1| ). In particular, for 

h s 1 — s 



< 1 



(4.3) 



and 



« fi? ^ ^ 



« ^1 ^2 



1 - g g 



h 



one gets: 



e « 1 



(4.4) 



(4.5) 



(4.6) 



The initial spin comes to its end when the (negative) quantity — 2 /2 Tkin approaches 
zero, and k ^ 1 . The next, second stage will be precession in the vicinity of separatrix 
(though without crossing it yet). Crossing of the separatrix may result in chaotic flipovers. 
The third stage is that of tumbling. It begins when (J^ — 2 I2 Tkin) = and k = 1 . 
It ends when (J^ — 2 I2 Tkm) > and k is smaller than unity, though not 
small enough to approximate the Jacobi functions by trigonometric functions. Mind that 
the transition from the initial spin to tumbling leaves parameters a , /3 , 7 , u and k 
continuous: none of these undergo a stepwise change. The fourth stage will be called 
"the final spin" : it takes place when the relaxation is almost over and A; ^ . At this 
stage (as well as during the initial spin) the Jacobi functions may be well approximated 
by trigonometric functions. 

As explained in the previous section, the suggested scenario is, of course, too ideahsed: 
on the one hand, small occasional interactions will easily force a nearly-prolate body to 
rotate around its minimal- inertia axis for most time; on the other hand, the relaxation 
rate will vanish in the vicinity of pole C , so that the perfect relaxation will never be 
achieved. 



V. DYNAMICS OF A FREELY PRECESSING BODY. RELAXATION RATE. 

We are interested in the rate at which a freely spinning body changes its orientation 
in space, i.e., in the rate of alignment of the maximal-inertia axis along the (conserved) 
angular momentum. 

In the case of an oblate body (Is > I2 = h), one could start with a trivial formula 
d6/dt = {dTkin/ dOY^ [dTkin/dt) where 6 stands for the angle between the major- 
inertia axis and the angular momentum (Efroimsky & Lazarian 2000)^^. This formula 
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would work since in the oblate case 6 remains almost unchanged through a precession 
cycle. Unfortunately, in the general case of a triaxial rotator, even in the absence of 
dissipation, this angle evolves in time. Luckily though, its evolution is periodic (formulae 
(Al) - (A4) in Appendix A), so that instead of using 9 one can use its average over a 
cycle. In practice, it turns out to be easier to operate with the time-average of sin^ 6 : 



d < sin 6 > d < sin 6 > dTi 



kin 



dt dTkin dt 

As shown in Appendix A, formula (AlO), 



(5.i: 



J2 <sm'^e> 



^2hT,,^ - J^) + h h Hi + ? + ^ + 0{k')\ . (5.2) 



h ^ ' h - h 1 \ 8 16 



A substitution of ( p.l4| ) for /? and k into the above formula entails (see Appendix A, 
equation (A12)): 

j2 ( d <sin^g> \ ^ h{hh + hh - 2/1/2) ^ 

+ 7 \Y'':~f { -A^ + oik^) . (5.3) 



4 J2 - 2/iTfe,„ V/3 - h) h - h 
while a substitution of ( ^.16| ) into ( |^ ) yields (Appendix A, formula (A15)): 
2 ( d <sin^g> \ ^ /1/3 _ 1 /i/3(/3 - h) - 2/iTfc,„ . 

\ dTkin Ji^nearA) ^3 " A 4 (/g - h) (h " A) 2 /a Tfc,„ - J2 ^ ^ 



Formulae ( [5.3| ) and ( |5.4| ) explain how the losses of the kinetic energy of rotation make 
< sin^ 9 > change. Since the kinetic energy decreases because of the inelastic dissipation. 



- kin 



W , (5.5) 



what we have to find is the rate of the elastic-energy losses W, quantity W being the time- 
dependent part of the elastic energy stored in the body due to the alternating stresses. 
Then, with aid of (|5?1|), ( |5.3| ), (|5.4|), we shall compute the rate of alignment: 



d < sin^ 9 > d < sin^ 9 > dW 
dt dTkin dt 

The next four sections will be devoted to the calculation of the dissipation rate dW/dt . 



VI. ESSENTIAL NONLINEARITY IN THE PRECESSION-CAUSED 
DISSIPATION. CASES OF HOT AND COLD BODIES. 

Our goal now is to describe the kinetic-energy dissipation caused by the deformations 
of the body, experienced in the course of its precession. The deformation of body is neither 
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purely elastic nor purely plastic, but is a superposition of the former and the latter. It 
is then to be described by the tensor e.^ of viscoelastic strains and by the velocity tensor 
consisting of the time- derivatives e^ . The stress tensor will consist of two components: 
the elastic stress and the plastic (viscous) stress. In the simpliest, so-called Maxwell- Voigt 
model, the components are additive (Tschoegl 1989)^^: 

= 4^ +4^ , (6.1) 

where the components of the elastic stress tensor are interconnected with those of the 
strain tensor (Landau and Lifshitz 1970)^^: 

eij = S,, + ( 4) - I 5,, Tr ^ , (6.2) 



a^f = K5ij Tre + 2 // (e,,- - ^ 5ij Tr , (6.3) 

/X and K being the adiabatic shear and bulk moduli, and Tr standing for the trace of a 
tensor. Components of the plastic stress are connected with the strain derivatives as 

^ij - Sij + {4^ - I Sij Tr ^ , (6.4) 



C 5ij Tre + 2ri (e,,- - ^ 5,,- Tr e) , (6.5) 



where rj and ( are the shear and stretch viscosities. 
Dissipation may be taking place at several modes: 

Q{uj) being the so-called quality factor of the material, and Wo(a;n) and < W{ujn) > 
being the maximal and the average (over a period) values of the appropriate-to-c<j„ fraction 
of elastic energy stored in the body. The average (over the precession cycle) of the total 
elastic energy reads 

<W>^ ^ J dV <aij eij > , (6.7) 
and it must be decomposed in a sum over the frequencies: 

<W>^Y. < > , (6.8) 

n 

For example, in the case of a symmetrical oblate body, studied in (Lazarian & Efroimsky 
1999) and (Efroimsky & Lazarian 2000), both the stress tensor a and the strain tensor e 
contain only the precession frequency a; . Therefore their contraction cTjj eij contains two 
frequencies: u and 2uj , and hence in this case dW/dV = dW'^^^ /dV + dW^'^'^^ /dV . 
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All in all, the general expression ( p.6|) entails: 




Un d < W{uJn) > 



uj„) dV 



\ 



(6.9) 



where the integral is taken over the entire volume V of the body. In the latter expres- 
sion we have deliberately put the quality factor under the integral, implying its possible 
coordinate-dependence. The coordinate dependence of attenuation should be taken into 
account whenever one is dealing with precession of an inhomogeneous body. We mean, 
for example, the problem of rotational stability of a spacecraft. Wobble of a strongly 
inhomogeneous asteroid is another example of relevance of the coordinate dependence of 



to be a higher-order correction. Instead, it is the higher-than-w frequencies 
that contribute the overwhelming share of the entire effect. This, crucial circum- 
stance had gone unnoticed in the preceding studies (Burns and Safronov 1973)^^, (Purcell 
1979), and was studied only this year in (Lazarian & Efroimsky 1999)^^ and (Efroimsky 
& Lazarian 2000)^^. In the latter two articles we were considering a simple case of a 
symmetrical oblate body (Ji = 12). In that case, the second mode was generated due to 
the quadratic dependence of the centripetal acceleration upon the angular velocity Vt : 
since the angular velocity of an oblate body precesses at a rate uo = [h — 1)VL^ (where 
h= /3//1 = J3/J2), the emergence of double-frequency terms in the expression for accel- 
eration (and therefore, in the expressions for stresses, strains, elastic energy, and finally, 
in the expression for the relaxation rate) is unavoidable. For the first time the presence 
and role of the double-frequency terms was discussed in (Lazarian & Efroimsky 1999)^^, 
in the context of cosmic-dust alignment, and in (Efroimsky & Lazarian 2000)^"^, in the 
context of cometary and asteroid wobbling. It turns out that this second mode often gives 
a leading input into the dissipation process. This is an example of a nonlinearity giving 
birth to a leading-order effect. It remains a puzzle as to why this, leading, effect had not 
been studied thitherto. It would be though unfair to say that the effect had gone com- 
pletely unnoticed. After our two articles had been published, Vladislav Sidorenko drew 
our attention to the fact that the second mode had been mentioned back in fifties by 
(Prendergast 1958)^^ and then forgotten. Prendergast took into account the centripetal 
acceleration but missed the term 17 x r in his analysis. He also ignored the emergence 
of the higher modes. Anyway, we would credit Prendergast for first noticing the nonlinear 
nature of the process. 

In 1973 Peale published an article (Peale 1973) devoted to inelastic relaxation of nearly 
spherical bodies, where he did take the second harmonic into account. 

Since in the case of an oblate symmetrical body only two modes are present, formula 



for W simplifies a lot: WoUate = W{uj) + W{2uj) = -2[uj < W{uj) > /Q{uj) + 



2u < W{2u) > /Q(2cu)] ^ -2ij < W{ij) + 2W{2uj)> /Q{to) where we used the fact 
that the quality factor depends upon the frequency very slowly: Q{2ui) ^ Qiyj) . This 
neglection of the frequency- dependence of Q is certainly valid when we consider inputs 
from frequencies differing from one another by a factor of two. However, in the generic 
case, when a broader band of frequencies comes into play, the frequency-dependence of 
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the quality factor in (|6.6| ) must be respected. This dependence may be crucial if the body 
is a composite structure with resonant eigenfrequencies of its own: attenuation at these 
may be especially effective. 

As we shall see below, whenever the rotation axis is (almost) parallel either to the 
maximal- or minimal-inertia axis of the body, the dissipation is taking place on two 
frequencies solely: uo and 2uj . This situation will be reminiscent of the above-mentioned 
case of a symmetrical oblate body. Relaxation in the vicinity of the separatrix is a far 
more complicated case; in that case numerous frequencies will be generated, and the 
frequency dependence Q{'jj) will be relevant. 

The quality factor Q is empirically introduced in acoustics and seismology to make 
up for our inability to describe the total effect of a whole variety of the attenuation 
mechanisms (Nowick and Berry 1972)3°, (Burns 1986)3^ (Burns 1977)^^ (Knopoff 1963)^1 
A discourse on the frequency- and temperature-dependence of the Q-factor is given in 
Appendix B. 

Another issue worth touching here is that of elasticity and plasticity demonstrated by 
materials at various temperatures. In order to calculate the terms entering ( |6.8|) , one must 
know the stress tensor (that can be found from knowing the acceleration of an arbitrary 
point of the body) and the strain tensor (that depends upon the stress tensor through the 
system of equations (|6.1D, ( |6.3D and (|6.5|)). In general, it is difficult to resolve the system 
( |6.1| ), ( |6.3|) and ( |6.5| ) with respect to eij{x, y,z). Fortunately, in two simple practical 
cases the system solves easily. These are the cases of cold and hot (plastic) body. 

As well known, at low temperatures materials are fragile: when the deformations 
exceed some critical threshold, the body will rather break than flow. At the same time, 
at these temperatures the materials are elastic, provided the deformations are beneath 
the said threshold: the sound absorption, for example, is almost exclusively due to the 
thermal conductivity rather than to the viscosity. These facts may be summarised like 
this: at low temperatures, the viscosity coefficient r] has, effectively, two values: one 
value - for small deformations (and this value is almost exactly zero); another value - 
for larger-than-threshold deformations (and that value is high. (Effectively it may be put 
infinity because, as explained above, the body will rather crack than demonstrate fluidity.) 
Therefore, within the range from the absolute zero up to at least several hundred degrees 
of Celsius the plastic part of the stress tensor may be well neglected. 

At high temperatures materials become plastic, which means that the shear viscosity 
rj gets its single value, deformation-independent in the flrst approximation. On the one 
hand, this value will be far from zero (so that the scattering of vibrations will now be 
predominantly due to the viscosity, not due to the thermal conductivity). On the other 
hand, this value will not be that high: a plastic body will rather yield than break. All this 
is certainly valid for the stretch viscosity ( as well. As a result, at temperatures higher 
than about three quarters of the melting temperature one may neglect the elastic part 
of the stress tensor, compared to its plastic part. (I am deeply thankful to Shun-Ichiro 
Karato for a consultation on this topic.) 

To simplify the stress tensor, we model the body by a rectangular prizm of dimensions 
2a X 2b X 2c. The tensor must obey three demands. First, it must satisfy the relation 

diCij = p aj (6.10) 
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aj being the time-dependent parts of the acceleration components, and paj being the 
time-dependent parts of the components of the force acting on a unit volume. Second, 
tensor aij must be symmetrical and, third, it should obey the boundary conditions, i.e., 
the product of the stress tensor and the normal unit vector, (Jijrij should vanish on the 
boundaries of the body (this condition was not fulfilled in (Purcell 1979)^^). 

It would be important to emphasise that the above assumption of the body being a 
prism brings almost no error into claculations performed for real irregular-shaped physical 
objects, like asteroids or cosmic-dust grains. The reason for this is that an overwhelming 
share of dissipation is anyway taking place not near the surface but in the depth of 
the body. This is especially evident from formulae (|8.71 ) - (|8.171 ), and it is just another 
manifestation of Saint- Venant's principle of elasticity. (I am grateful to Mark Levi for 
drawing my attention to this fact.) So, whether the body is indeed a rectangular prism 
or more like an ellipsoid, will not make much difference for an estimate of the relaxation 
time. Mind though that for shells Saint- Venant's principle does not work, so that in the 
case of spinning spacecrafts the subtleties of their shape may be relevant. 



VII. THE ELASTIC ENERGY OF ALTERNATE DEFORMATIONS 

Unless the temperature is too high, the bodies manifest, for small deformations, no 
viscosity ( ujt] ~ ^ fi ^ K ), so that the stress tensor is approximated to a very 



high accuracy by its elastic part: instead of the system (|6.1| ) - ( |6.5|) one may write 



en = 6ii + ( an — - da Tr (x) , (7.1) 

This will enable us to derive an expression for the elastic energy stored in a unit volume 
of the precessing body: 

d <W > 1 1 f f2fi 



dV 2 ' ' 4/i [\9K 



- < (Tr > + < aij aij ^ 

I - Y^—I < > + <(^L> + <(^yy> + < ^Iz > < C^ly + CTyz + ^^L 



1 r 1 
4/i 



where we have made use of the expressions connecting the shear and bulk moduli with the 
Young modulus and Poisson's ratio z/: since K = _E/[3(1 — 2z/)] and fi = E/[2{l + h')] 
then 2fj,/{9K) — 1/3 = — i//(l + z/) . As Poisson's ratio u is, for cold sohds, typically 
about 0.25, one may safely put 2fi/{9K) — 1/3 ~ — 1/5 . 

The above expression (|7.2| ) must be decomposed into a sum like (|6.8| ) . In what follows 
we shall be interested in the energies averaged over the precession period. For this reason, 
in (|7.2| ) all af, 's are averaged. (See Appendices C and D for details.) 
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VIII. RELAXATION RATE IN THE VICINITY OF POLE C: 
THE RELAXATION IS ALMOST COMPLETED AND THE BODY IS 
SPINNING ALMOST ABOUT ITS MAXIMAL-INERTIA AXIS 



Near the poles parameter k is close to zero. This justifies the following simple asymp- 
totics for elliptic functions (Abramovitz & Stegun 1964, formulae 16.13.1-3)^^: 

cn{u, k'^) = cos u + — k"^ [u — sinu cosu) sinu + O(fc^) , (8-1) 



4 

sn{u, k"^) = sin u — -k'^ {u — sinti cosm) cosm + 0{k'^ 
dn{u, k^) = 1 - sin^M + 0{k^) , 



(8.2) 



13) 



In the vicinity of pole C, i.e., during the "final spin" (when fl is almost aligned along 
or opposite J ), we substitute ( |2.12[ - p^.l^ and ( |2.19| - |2.27| ) into ( ^.181 ). Then we must use 
the asymptotics ( p.l|) - ( ^.3|) neglecting terms of order higher than k"^ . Mind that, for 
J2 > 2Tkinl2 , the parameters P and 7 are of same order as k (as evident from 
( p.l4| )), while R is of order k^ (according to (pl8|)). This will give us: 



.{*) 

(near C) 



ei 



{ - X [P + (1 



Q) 0^ sn^[u, k"^] + y (^a uj k"^ + P sn[u, k'^] cn[u, k"^] + 
+ z (Puj + ay) cn[u, k^] dn[u, A;^]| + 



+62 [x (^py - aujk'^^ sn[u, fc^] cn[M, A;^] - y P + R - (Q + S) p'^ sn'^[u,k'^ 



+ 



+ z {a P + 70;) sn[M, /c^] dn[u,k^]^ + 

+63 |x( — /So; + ay) cn[u, /c^] dn[u, k"^] + y (a P — y 00) sn[u, k"^] dn[u, fc^] — 

- z [r+ [l - S)P^ sn2[u, k^]] } = (8.4) 

61 I - X (1 — Q) cos 2m + - ?/ {a uj k'^ + P y^ sin 2m + z [P uj + ay) coswj + 
+62 |-x {Py — acuk'^^ sin 2m — y -{Q + S)P'^ cos 2m + z{aP + yu) siuMj + 
+ 63 ^x { — P u + ay) cos u + y [a P — y uj) sin w + z -{\ — S) P"^ cos 2m | + 
+ I time — dependent k"^ — order terms originating from the 
k^ — order terms in expansions (8.1) — (8.3) | + 



+ { time — independent terms } 
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where u = uj {t — to) . In the further calculations we shall ignore the time-independent 
terms emerging in ( p.5[ ) because, in order to calculate the inelastic-dissipation rate, we 
need only time-dependent part of the stress tensor. 

Dissipation is taking place in two modes one of which has the frequency of precession, 
while another one is of twice that frequency. (If one plugs into (|8.4|) all the high-order 
terms from ( |8.1| ) - ( ^.3[ ) they will give an infinite amount of the higher modes in ( |8.4| ). In 



the vicinity of poles, we neglect the high-order terms in ( |8.1|) - ( p.3| ), and thereby neglect 
harmonics higher than second.) As already mentioned in Section VI, the second mode 
originates from the centripetal term in ( p.l7| ). This fact is understood especially easily if 
we assume that the body is oblate and symmetric. In this case one component of fl (the 
one parallel to the axis of maximal inertia) will stay unchanged, while the other two will 
be proportional to sin ut and cos ut (which simply means that f2 is precessing at rate 
u about the maximal-inertia axis). Quite evidently, squaring of f2 in ( p.l7|) yields the 
double frequency. Mathematically speaking, in the case of an oblate body the realm of 
applicability of the solution ( p.l5| ) - ( p.l6| ) shrinks to a line, so that the solution ( |2.13| ) 



( p.l4|) accounts for the entire process. As evident from ( p.l4| ), in the oblate case A; = 



and therefore formulae ( ^.11 ) - ( |8.3| ) contain only terms of order . 

So we shall strip ( |8.4| - p.5|) off its time-independent terms, and shall plug the 
time — dependent terms into ( [6.101 ). Integration thereof will then give us expressions 



for the time — dependent components of the stress tensor: 

= - I (1 - Q) {x^ - a^) {sn^[u;{t - to), k^] - < sn^[u;{t - to), k^] > } = (8.6) 



2 ^ 



Q) - a') cos2[u;(t - to)] + 0(A;2)| , (8.7) 



'^yy 



^ {S + Q) 13^ (y^ - b^) {sn^[uj{t - to), k^] - < sn^[uj{t - to), k^] > } 



2 



^(S + Q) {y' - b') {- ^ cos2[c.(t - to)] + 0(P)} , (8.9) 



0", 



^ (1 - ^) (z^ - c^) {sn^[u{t - to), - < sn^[u{t - to). A;'] > } = (8.10) 
= -^{1- S)P' (z^ - c') |-lcos2[.;(t - to)] + ^(fc^)} , (8.11) 



^ [{p^ + aujk^) {y'-b-") + 
+ (/57 - atufc^) (x^ - a^)} [sn[uj{t - to), k"^] cn[uj{{t - to), k'^] - 



< 



sn[uj{t — to), A;2] cn[u;((t-to), A;^] > } = (8.12) 
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+ (/37 - aue) [x^ - a^)} |i sin2[cu(t - to)] + 0(A;2)| , (8.13) 



P 



+ {-/3(jj + a-f) (x^ - a^)} {dn[uj{t - to), k'^] cn[uj{t - to), A;^] - 

— < dn[oo{t — to), k ] cn[uj{t - to), > } = (8.14) 

= ^{{(3u + a^) {z'-c') + 
+ {-puj + a-f) (x^ - a^)} [cos[uj{t - to)] + 0(P)} , (8.15) 

'^yz = I {{ap + uj^) [z^ - c^) + 
+ {a/3 - CU7) {y^ - h^) } [dn[uj{t - to), k"^] sn[uj{t - to), k"^] - 

— < dn[uj{t — to), fc^] sn[cj(t - to), fc^] >} = (8.16) 

= P-{{aP + u^) {z'-c^) + 

+ (a/3 - u^) (y^ - h^)] {sin[a;(t - to)] + 0{k^)] , (8.17) 

The symbol < ... > stands for averaging over the mutual period r of functions sn 
and cn . (See Appendix A.) In the above expressions, it might be better to write afj 
instead of cjij , in order to stress that we are considering only the time-dependent part, 
but we would rather omit the superscrips for brevity. The above expressions ( ^.5[ ), (|8.7| ), 
(IJ), ( CT ), (|87[3| ), (|87[5D , ( gATl) coincide, in the limit of oblate symmetry {h = h). 



with formulae (19) - (23) from our previous article (Lazarian and Efroimsky 1999)^^. 

The expression (8.5) is exact, while the formulae (8.6) - (8.17) implement the poly- 
nomial approximation to the stress tensor. This approximation keeps the symmetry and 
obeys (6.10). The boundary conditions are satisfied only approximately, for the off- 
diagonal components. This approximation very considerably simplifies calculations and 
yields only minor errors in the numerical factors (8.26) - (8.28). A comprehensive analysis 
of the polynomial approximation will be presented elsewhere. 
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To calculate the dissipation rate, we shall need averaged over the precession period 
squares of the above stresses, < afj > , as well as < {Tra)'^ > . Moreover, for our goals 
we shall need these calculated up to terms of order k"^ inclusively. This demand makes 
it necessary to have the above expressions (g^), (^), (pl| ), (|87[3| ), (|87[5| ), (|8T? 



with all the fc^— order terms written explicitly. How to get these terms? On the face of 
it, the answer is trivial and looks like this. In the above formulae we approximated the 
elliptic functions using only order terms of 18.3|) ; now, let us keep also the /c^-order 
terms. Surprisingly, this is the case when the simpliest shortcut leads to a wrong answer. 
Plugging of the contained in (U - U) fc2_order terms into (|8]3), (|]|), (|8lll) . 



!T2l) , ( ^.14| ), (|8.16|) , with the further squaring thereof, will give birth to secular terms in 
the expressions for < afj > , i.e., to terms linear in m = uj(t — to) . Averaging of these 
terms will entail ambiguities: one will get into an illusion that it does matter whether 
to integrate from through 2 vr or, say, from 2 vr through 4 vr . The secular terms 
have been long known in nonlinear mechanics and astronomy where they often tarnish 
calculations and sometimes become a real pain. Luckily, in our case we can sidestep this 
obstacle by employing directly the fundamental definition of the elliptic functions: 

1 /2 

sn(u, k^) = sin0 , cn(M, k"^) = cos(f) , dn(u, k"^) = (^1 — fc^ sin^ 0^ (8.18) 
the auxiliary quantity being connected to u like that: 

u ^ / . (8.19) 

^0 fi - p ^i^^eV 



This will give us the key to a correct calculation of the averaged-over-period quadratic 
and quartic forms. For example, the average < sn^w dn^u > will read: 

< sn^iu, k'^) dn^iu, k"^) > = — f sn^iu, k"^) dn^iu, k"^) du 



r Jo 



• 2 , /i ,2 ■ 2A du ,^ 1 /■^'^ . 2 , A ,2-2 ,\l/2 

sm (p [1 — k sm m —— dm = — — / sm (p [1 — k sm 




AK 27r Jo V 2 W 1 + k^/A ^ 2 



where we used (A7). The squared and averaged in the above manner stress components are 
presented in Appendix C, expressions (C2), (C4), (C6), (C8), (CIO), (C12) and (C14). 
Substitution thereof into ( [7.2|) will lead us to the expression for dissipation per unit 
volume: 

d <W> d <iy(2.^)> d <iy('^)> 

dv - dv + dv ^^-^^^ 

where the first term stands for the dissipation of oscillations at frequency 2 u : 

d < 1^(2<^) > 
Jt 

{ - I +Vi < ^^^^)' > + < > + < > + < > + 2 < < > } , (8.22) 
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while the second term expresses the dissipation at the principal frequency: 



dt = 4^ {2 (<<> + <-i>)} (8-23) 

Expressions for d < W^'^^ > / dV and d < W^'^^^ > / dV in terms of /i,2,3 are presented 
in the Appendix (formulae (C20) and (C21)). These expressions should be now multiplied 
by 2uj/Q{(jj) and 4:Uj/Q{2uj) , correspondingly, and integrated over the volume of the 
body, as in ( |6.9|) . The outcome of this integration will be the total dissipation rate W 
that must be plugged, together with ( p.3[ ) into (|5.1|). Here follows the result: 



d <sin2^> 2/3p2 {2hTkin - 3' 



dt /i Q{uj) J2 

+ 



LU (j^ — 2IiTkin) Hi 



ih - h) ih - h) 



1 h{h-hf 2/3Tfc,„-j2 



4 {h-h) [h-hf J' - 2/lT; 



kin 



UJ Ho [2 I3 Tkin — 



The ratio Q{ijj) /Q{2uj) is typically close to unity, unless the structure of the body or 
the properties of the material provide resonances. Terms with Hq and Hi are due to 
the dissipation of oscillations at frequency uj , while the term with H2 is due to the 
vibrations at 2uj . 

Numerical coefficients Hq , Hi and H2 emerging in ( |8.24| ) are geometrical factors 
that depend upon the moments of inertia and dimensions of the body. General expressions 
for Hq 12 are given in Appendix C. Obviously, Hq vanishes in the oblate case. 

Equation ( ^.241 ), together with ( ^.141 ) and with equation 

J2 < sin^^ >(_c) = \ ^ tvr"^!'/' (2^3T..„ - J^) + 0{k') (8.25) 

connecting Tkin with < sin^ 6 > , makes a system of equations describing relaxation 
in the vicinity of pole C. (Equation (|8.25|) is a truncated version of (|5.2| ). For details see 
(All) in Appendix A.) 

Let us elaborate on the factors Hq , Hi and H2 ■ In the case of a homogeneous body 
of dimensions 2a x 26 x 2c , expressions for the factors read (see the end of Appendix 
C): 

^ ^ 317 ahl / . a' \ .fil 

' m4 (62 + c2) (a4 - c4) {a^ + b^) [b'' - ^ - d J ' ^ ' 



H2 : 
and 



100 a^b^c - aH^c^ - 
m4 Ic 



'b^e> + 0.21 6 + 0.19 6^ + Q.21 ab'^c^ 



+ b^y 



(a4 - c^Y (64 _ c4) 



(8.27) 
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The denominators of (C22) - (C24) contain expressions {I3 — I\ 
result, the denominators of ( ^.26| ), ( p.27| ) and ( |8.28| ) contain [a^ — c 
It would be appropriate to make sure that nothing wrong happens when a 



and (/a — J2) ; as a 
and {h" - d'f . 



c or 



b . We assumed from the beginning that I3 > h ^ h , i-e., that (for a prism) 



a > b > c. Therefore it would be enough to investigate the case of a 
that the parameter k given by (|2.14|) and ( p.l6|) never exceeds unity: k 



and = 1 at the separatrix. From ( |2.14| ) we see that, for {I^ — I2 



^ b . Recall also 
= at the poles 
(62-c2) ^ 0, 

the condition k < 1 is fulfilled only if fi^/^l < (I^^I^ - I^)) / {h{l2 - h)) . In other 
words, making {I3 — I2) (and (b — c) ) infinitesimally small leads to infinitesimal squeezing 
of the region around pole C between the separatrices on Fig.l. Thus the region, where 
the appropriate solution is applicable, shrinks into a point. 

In our analysis it is possible to get rid of the variable Tkin completely: one should 
express it through < sin^6' > by means of ( |8.25| ), and plug the result into ( |8.24| ). This 
will give us what we would call Master Equation, a differential equation for < sin^ 9 > : 



d < sm 9 > 
Jt 

1 - 2 < sin^ 9 > 
1 

+ - 



h {I3 - h) 



< sin^ 9 > 



hh + hh - 2hh_ 

ih - hfll <sin2^> 



[/1/3 + hh - 2hh 



2 (/1/3 + hh - ^hh) - 2 < sm'9> h ih - h 



LU Ho2h ih - h) < sin^ 9> + 



+ 4u;H2^^ <^n'9> hih - I 
Qiu) 




where, according to ( p.l4|) and (|8.25|) 



'h 



\ 



{h-h)ih~h 
hh 



2 < sin^ 9 > 



h ih - h 



hh + hh - 2/1/2 



^.29) 



^.30) 



Equation (|8.29|) is one of the main results of our study. It describes the relaxation in the 
vicinity of pole C corresponding to rotation about the maximal-inertia axis. Simply from 
looking at this equation one can understand several important features of the relaxation 
process. To start with, it follows from ( fj.29| ) that d < sin^ ^ > /dt vanishes in the 
limit of (/a — /i) ^ , which naturally illustrates the absence of relaxation in the 
case of all moments of inertia being equal to one another. Second, the overall factor 
< sin^ 9 > standing before the brackets in the right-hand side of ( |8.29|) evidences of a 
gradual decrease in the relaxation rate: the major-inertia axis will be approaching the 
angular momentum vector but will never align along it exactly. 
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Technically, the Master Equation (|8.29| ) becomes a self-consistent differential equation, 
describing the time-evolution of < sin^ 9 > , only after the expression ( |8.3CI| ) for uj is 
plugged into it. We did not bother to do this not only for the sake of brevity. In fact, 
equation ( p.29| ) as it stands is of more practical interest than the self-consistent differential 
equation. It enables, for example, an astronomer to use the measurable quantities u and 

< sin^ 6' > , to predict the relaxation rate in the short run. In the real life "short run" 
means: the time span during which the currently available resolution of the optical or radio 
equipment makes it possible to notice the narrowing of the precession cone. Nowadays 
spacecraft-based equipment provides an angular precision of 0.01° and even better. This 
gives us a chance of observing precession damping within a period varying from several 
months to several years, for different objects (Efroimsky 2000)^^. Soon Rosetta mission 
will give it the first try (Hubert & Schwehm)^^, (Thomas et al)^^. 

Now let us briefly dwell on the limit of an oblate body { I2 = h )■ In this case, 
the precession is known to be circular (Efroimsky & Lazarian 2000)^^, so the averaging 

< ... > may be omitted. The simplified Master Equation will then look: 



d (sin 



(oblate) 



dt 



43' h-h .2 
sm 



+ 2uH2 



e [ujHih cos2^+ 



(8.31) 



where 



(jj = — I — — 1 1 cos (J . 

h V-'i 



(8.32) 



For an oblate homogeneous rectangular prism, the latter and the former, with ( ^.26|) and 
( ^.27| ) plugged in, will give: 



3 • 3^ 

- sm^^ 



63 ( - cot^ ^ + 20 
aj 1 



l-2{c/a)^ + 0M{c/a) 



2 (c/a)4 + (c/a)8 



/■ \ (oblate) 
Jt) 



pQ{uj) 



[1 + (c/af 



— - sm'^6' 
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63 {c/af cot^e + 20 



[1 + {c/a 



214 



where 



P Q{oj 



^ + 0{k\ {c/af) (8.33) 



(8.34) 



and it is assumed that Q{uj) ~ Q{2uj) . This perfectly coincides with the exact formula 
obtained in (Efroimsky & Lazarian 2000) by a rigorous treatment possible in the oblate 



case: 
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.^.ioUate) ^ ^^^^^ 

\ / {EkL 2000) 



63 {c/aY cot2 + 20 



[1 + {c/a 



214 



, 8.35 



Now let us see what happens with the Master Equation ( ^.291 ) when the shape of the 
body is almost prolate { I3 ~ I2): 



\ dt J 

fj'Qiuj) I V 2 

-LjHo2{h - I2) + ^H2u; < sin' 9> {h - h) ] + 0{k^) (8.36) 

Even though the term containing H2 contains also multiplier (Is — I2) , it diverges 
in the limit of prolate symmetry (see (C23)). However, there is nothing bad about it, 
as explained in Section III: one should not make (Is — I2) approach zero for fixed 9 , 
but rather fix some value of (Is — I2) , small but finite, and then make 9 decrease 
to zero. As already mentioned (see the comment after ( ^.27|) ), as the shape approaches 
the prolate symmetry, the applicability region of the solution shrinks. Still, the fact is 
that within the applicability region (called in Section III the final spin) the second-mode 
term will not necessarily be much less than the first one. The ratio of these terms will 
depend upon sin' 9 , which means that the typical time of relaxation may be a steep 
function of the angle. This typical time must be proportional, for dimensional reasons, to 
/i Q/(2 J' p' (/a — h)) , but the numerical factor may be quite 6*— dependent, due to the 
presence of the second term in ( ^.361) . We had to dwell on this subtlety due to its practical 
relevance. In the recent literature they sometimes use the formula for relaxation time, 
derived for oblate bodies, in order to estimate relaxation of a tumbling prolate rotator. 
We did this in (Efroimsky & Lazarian 2000)'^ when discussing asteroid 4179 Toutatis, 
while (Black et al 1999)'^ employed this estimation for asteroid 433 Eros. We see now 
that this was wrong even for the final spin about pole C. The more so, it was absolutely 
unjustified to use this estimate for a tumbling body (i.e., in the vicinity of the separatrix) 
as done in the said articles. 



IX. RELAXATION RATE IN THE VICINITY OF POLE A: 
THE BODY IS SPINNING ALMOST ABOUT ITS MINIMAL-INERTIA AXIS 



In the vicinity of pole A, i.e., during the "initial spin" (when the angular velocity 
f2 is almost perpendicular to the maximal-inertia axis, and ^ ), we substitute 
asymptotics d^jp) into (|J3|), (|J5| - [2A6| ), (|2l9| - |;23D and - |3T]), the results 
to be plugged into (|2.18|). This leads to the following expression for the acceleration: 



{near A) 



X 



P + {l-Q)P 



,2 2 
sn u 



+ y (a u; + /? 7) 



sn[u, k 



dnlu, k'^ 
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+ z {(3uj + a'f) cn[u, k^] dn[u, k"^] | + 



+ 62 {x - auj) sn[u, dn[u, k'^] - y P + R - {Q + S) (3"^ sn^[u, + 

+ z {a (3 + '-) uok"^) sn[u, k"^] cn [u, k^] | + 

+ 63 |x ( —(3 uj + a'-)) dnu cn[u, k'^] + y {a (3 — 70; fc^) sn[u, k^] cn[u, k^] — 

- z[r +(1 - S)(3^ sn^[u, k^]]} = (9.1) 

= ei I - X (1 — Q) cos 2u + y {a uj + /5 7) sin u + z {j3 uj + a '-/) cos m| + 

+ e2 I X (/5 7 — a ) sin -u — - {Q + S) y f3'^ cos 2u + -z {a f3 + 'y u k'^) sin 2m| + 

+ 63 l^x { — l3uj + a'-f) cos u + -y {a (3 — 'J UJ sin 2m — z - [1 — S) /3'^ cos 2m| + 

+ I time — dependent terms of order k'^ ^ + { time — independent terms } . 



En route from (|9.1| ) to (|9.2| ) we employed asymptotics (|8.1|) - ( p.3| ), then separated out 
the time-independent terms (which may be dropped, because they do not influence the 
inner dissipation), and we also neglected terms of order higher than fc^ (we remind that, 
according to (|2.16| ), for < 2Tfcj„ J2 , the parameters /3 and a are of order k). In 
the above expression, u = oj (t — t^) . Parameters S and Q are expressed by (|2.8D. 
Parameters a , /3 , 7 , uj and k are expressed by ( |2.16| ) and thus are different from a , 
/3 , 7 , and k used in the preceding section (where they were expressed by ( ^.14| )). 

Similarly to the preceding section, we shall use equation ( |6.10| ) and expression (|9l^ ) 
to compute the stress tensor. This will lead us to: 



^(1 



Q) f3^ (x^ - a^) {sn^[uj{{t - to), 



- < sn^[uj{{t - to), ^1 > } 
i p (1 - Q) (x^ -a') I - i cos2[^(t - to)] + ^(fc^)} , 



(9.3) 
(9.4) 



P 
2 



{S + Q)(3' {y' -b') {sn'[u{{t-to),e] 
- < sn^[Lu{{t - to), A;2] > } 
^p(S + Q)(3' (y'-b') 1^-^ cos2[^(t - to)] + 0(fc2)| , 



(9.5) 
(9.6) 
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= - ^ (1 - S) (3' {z' - c') {sn'[u;{{t - to), k'] - 

- < sn^[uj{{t - to), fc^] > } = (9.7) 
^ - lp(l- S)P' {z' -c') [-^ cos2[c.(t - to)] + 0(^2)| , (9.8) 



+ (/37 - auj) (x^ - a^)} {sn[tt;((t - to), k'^] dn[uj{{t - to), k"^] - 

- < sn[u;{{t - to), k'^] dn[u;{{t - to), k"^] >} = (9.9) 

= ^ P {(/57 - a^) - «') + (/37 + a^) (l/' - {sin[a;(t - to)] + 0{k^)} , (9.10) 

(Jxz = I {{P^ + «7) - c') + 
+ {- Pu + a-f) (x"^ - a^)} {dn[a;((t - to), k'^] cn[p((t - to), /c^] - 

- < dn[a;((t - to), k"^] cn[p((t - to), k"^] >} = (9.11) 
= ^ p {(/?a; + a7) (-2' - c') + (-/3a; + aj) (x^ - a^)} {cos[u{t - to)] + 0{k^)} , (9.12) 

= f {(a/3 + c^7fc^)(.^-c^) + 
+ (a(3 - uj-fk"^^ {y^ - 6^)} {cn[a;((t - to), k'^] sn[a;((t - to), A;^] - 

— < cn[Lj{{t — to), ] sn[a;((t-to), k^] >} = (9.13) 

+ (a(3 - uj^k') (y^ - h^)} sin2[c.(t - to)] + 0(fc')} , (9.14) 
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wherefrom we obtain the (averaged over the precession period) quantities that emerge in 
the expression ( [7 .21 ) for the dissipation rate: af^ and {Tr a)^ . These expressions are 
written down in Appendix D. Substitution thereof into ( [7.2|) , with the further integration 
gives the total energy W of alternating stresses. 

Similarly to (8.6) - (8.17), the above formulae give a polynomial approximation to the 
stress tensor. (See the comment after formula (8.17).) 

All the further scheme of calculation exactly repeats that from the preceding section. 
Energy W consists of two components, one on the principal frequency, another on the 
second mode. These should be multiplied by 2u/Q{uj) and 4 u;/Q(2co') , correspondingly 
(as in formula ( |6.6D ). It will give us the overall dissipation rate W . Plugging this rate, 
along with (p.4|) into (|5.6|) yields: 



d < sin 6 > 



T, = - T2 - 2 h Tu,n) {UJ (2h Tkin - J') 5i [1 

at fi Q[ijj) I3 — Ii ^ / L V / 



kin 



4/2-/12/3 T, 



kin 



UJ J2 - 2/1 Tkrn] + 2UJ 



Q{2uo) 



2 /l Ti^in S2 



where the geometrical factors 5*0 , Si and S2 are given in Appendix D. 

In the case of a homogeneous rectangular prism of dimensions 2a x 26 x 2c 
factors ^0 , Si and 5*2 read (see Appendix D): 



the 



(9.15) 



Si = 



86.5 



ah c 



(a2 + 62) (62 + c2) 



0^ + 1.70^4 + 68 a"" + 1.7 a^ + 
- + — 



a^ - 64 



a^ - c4 



, (9.16) 



^ _ 21 a 6c 

' " ^ (62 + c2)' (a4 - c^)' (a^ - h^f 



aH^ + 2.8 aH^c^ - A.^aH^c" + a«c^ 
-A.SaH''^ + 4.8 6«c« 



and 



5n 



32.4 



ah c 



+ 1.67 a^c^ + c8) (62 - c2) 



64) (a4 



C4)2 (62 + c2) 



(9.18) 



(9.17) 



As explained in the end of the preceding section, multipliers like {I^ — Ii) and {I2 — /i) 
in the denominators of the expressions for 5*0, Si , S2 presented in Apendix D, as well as 



A\2 



and 



6^) in the denominators of ( |9.16| ), ( p.l7| ) and ( |9.18| 



multipliers (a 
are harmless. 

Similarly to pole C, in ( |9.15D we have two contributions: the one with Sq and 5*1 
originates from the dissipation of oscillations at frequency uj , while the one with 5*2 
comes from 2uj . Often it is the second mode that dominates the dissipation. For the 
case of an oblate body {I^ > Ii = I2) this fact was proven in (Efroimsky and Lazarian 
2000)24. In the case of an almost prolate rotator, the importance of the second mode can 
be easily understood simply from looking at Fig. 1. We see that the trajectories described 
by the vector f2 remain more or less circular up to a close vicinity of the separatrix, 
i.e., that the trigonometric approximation of the Jacobi functions is valid through a fairly 
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large region. In this region, therefore, our formahsm does work. Let us estimate the input 
of the second mode at points D and F on Fig. 1. Point D depicts the situation when 
= and ^^2/^1 = 1/2 5 while point F corresponds to ^2/^1 = '^H ■ Following 
(Black et al 1999)^^ we have prepared the picture so that it corresponds to an example 
from real life, asteroid (433) Eros. To that end we assumed /2 = 3/i , h = 3.05 Ji , 
which is the same as a = 2.19 c , h = 1.05 c. A simple calculation using ( |9.16| - 19.17] ) 
and ( |2.5| - p.6|) shows that at point D the second-mode term in ( p.l5| ) is less than one 
tenth of the principal-mode term with 5*1 {So being negligibly small): 



2S2 

Si 



kin 



2 /, n 



0.08 



'3 kin " / £) 

At point F though, the second-mode contribution slightly dominates: 

2 5*2 / — 2 Ji Tfc, 



(9.19) 



- kin 



Si 



2 I3 Tkin — J" 



1.02 



(9.20) 



though in reality the nonlinear input at F is higher; first, because of the less-than-unity 
multiplier accompanying uj Si in ( p.l5| ) and, second, because of the higher-than-second 
harmonics. It is perhaps pointless to approach the separatrix closer than F, because 
the higher- frequency terms omitted in (|8.1| - |8.3|) will become relevant. Anyway, their 
relevance will only add to the nonlinearity. 

Equation ( |9.15|) , along with ( p.l6|) and equation 



y <sin'^> 



(near A) 



2 ^3 Fkin J 



2 /, 



y - 2hTk^n] + Oik"") (9.21) 



that follows from (|5.2|) (see also formula (A14) in Appendix A), constitute a self-consistent 
system of equations describing relaxation in the vicinity of pole A. By means of the latter 
equation, one may express Tkin through < sin^ 6^ > and plug the result into (|9.15|) . 
This would yield a differential Master Equation for < sin^ 9 > : 



d < sin 6 > 



VJ2 

// Q{uj) 



{I3 - h) (1- < sin^^ >) [ujSi 
1/3-/2 /i 



LO So 



2/1 



1 - < sin^ 9>] + 2ujS- 



2/2-/1 /3 

Q{2uj) 2/1 



dt 

2 < sin^ e> - 
^ ( 1 - < sin^ 9 > 



1 - 



Q(a;) I3 



1 - < sin 9 > 



(9.22) 



where, according to ( p.l6| ) and ( |9.21 



CO 




2 < sin^ 9> - I 



(9.23) 
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For sin^ 6* = , i.e., for f2 parallel to axis (1), the relaxation rate d < sin^ 9 > /dt 
vanishes. This means that vector will be leaving the unstable-equilibrium point (pole 
A on Fig.l) infinitesimally slowly. 

When (J3 — h) — > , the expressions (D18) - (D20) for S'0,1,2 diverge as (/a — 
. Therefore the right-hand side of ( |9.22|) will diverge as (Is — , instead of 

approaching zero as one might expect on physical grounds. The reason for this would- 
be divergence is that our analysis is applicable, as explained in Section II, only in the 
adiabatic approximation. Therefore ( |9.22| ), as well as ( ^.29[ ), works for as long as inequality 
(13) holds. 

Similarly to the Master Equation (|8.29|) , the Master Equation ( |9.22D becomes a self- 
consistent differential equation for < sin^ 6 > only after ( p.23|) is substituted into it. 
Similarly to (|8.29| ), here we have deliberately abstended from plugging ( |9.23|) into (|9.22|) , 
because without this substitution ( |9.22|) is of more practical use. (See the explanation in 
the end of the preceding section, between equations (|8.3CI|) and (??)). 



X. CONCLUSIONS AND PRACTICAL APPLICATIONS. 



Formulae ( p.22| ), for pole A, and ( ^.29|) , for pole C, constitute the main result of this 



article. These are differential equations for the relaxation rate of a precessing homogeneous 
body of arbitrary moments of inertia. The relaxation rate is defined as the rate at which 
the major-inertia axis approaches the angular momentum about which it is precessing. 

Formula ( |8.29| ) describes the relaxation rate of a body spinning approximately about 
its maximal-inertia axis (pole C), while (|9.22|) describes the relaxation of a body rotating 
almost about its minimal-inertia axis (pole A). 

In these formulae the contributions of the modes u and 2 uj are manifestly separated 
{u being the precession rate). Often the dissipation at 2u; gives a considerable input (as 
shown in the example ( |9.2CI| ) in the end of section IX), or even dominates, as in the case 
of oblate body, when Ii = I2 - see (Efroimsky & Lazarian 2000)^^*. 

Our formulae ( p.22| ) and ( |8.29|) were derived in assumption of parameter k being 



small. When does this assumption work? To get a simple answer to this question, let us 
look again at Fig. 1 . The approximation is valid on the part of the ellipsoid surface, 
that is covered with almost circular trajectories; a divergence of trajectory shape from a 
circle signals of inacceptability of the approximation (^[1]- |8.3| ). We see from the picture 
that, for example, for an almost prolate body our approximation remains valid not only 
in the vicinity of pole A but almost all way up to the separatrix. However, after the 
separatrix is crossed and the body begins librations, our approximation will regain its 
validity only in the closemost vicinity of pole C. Our formula (^.29) coincides, in the limit 



of oblate symmetry, with the exact formula ( |8.35| ) derived for oblate bodies in (Efroimsky 
& Lazarian 2000) 

The developed formalism has two immediate applications we know of. First of all, 
it is the study of wobbling asteroids (Harris 1998)^^. Wobbling may provide valuable 
information on the composition and structure of asteroids and on their recent history 
of external impacts. Astronomers observing precessing asteroids often ask: "Why do we 
have so few excited rotators in the Solar System?" (Steven Ostro, private communication). 
One of the reasons for this deficiency is that, due to the dissipation in the second mode 
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(Efroimsky and Lazarian 2000)'^'^ and higher modes, the effectiveness of the inelastic 
relaxation turns to be much higher that thought previously. 

Another obvious application of our formalism is the problem of cosmic-dust alignment. 
When they talk in astrophysics about the cosmic-dust alignment, they imply not the 
alignment of fl or the major-inerta axis along the angular momentum J but the 
alignment of the major-inerta axis relative to the interstellar magnetic field. This is 
a well-known phenomenon that manifests itself through the observable polarisation of 
starlight. There exists a whole bunch of physical mechanisms that make the cosmic dust 
align with respect to the magnetic field. Different types of these mechanisms dominate 
under different physical conditions, but all of them are based on the fact that cosmic-dust 
granules are swiftly precessing about the magnetic field lines. This precession is called 
into being by the interaction of the magnetic moment of the granule with the field. The 
magnetic moment is created by the Barnett effect and is thus parallel to the angular 
velocity f2 . Some of the known mechanisms of cosmic-dust alignment are very sensitive 
to the coupling between 17 and J , and this is when the inelastic relaxation comes 
into play. Until recently the Barnett relaxation was believed to be the leading relaxation 
mechanism. This viewpoint was expressed in the long-standing article by (Purcell 1979)^^. 
Remarkably, Purcell underestimated the effectiveness of the inelastic relaxation in the 
same manner as (Burns and Safronov 1973))^^ did it for asteroids: he missed the input 
provided by the second mode. Besides, he failed to satisfy the boundary conditions on the 
stresses. As a result, he underevaluated the effectiveness of the inelastic dissipation by 
several orders (see (Lazarian and Efroimsky 2000))^^). A study of the role of the inelastic 
dissipation in various mechanisms of cosmic-dust alignment is now on the way, and some 
results have already been published by (Lazarian and Draine 1999))^^. 

There may be a possibility that the developed formalism finds its third application in 
the research on spacecraft-rotation damping. 

XI. THE GENERIC CASE 

In the article thus far we have studied the dissipation in the vicinity of poles, i.e., 
the case when the body rotates about an axis that is close either to that of minimal or 
maximal inertia. All our formulae were derived up to 0{k'^) , parameter k being small 
near poles A and C (Fig. 1). 

In the general case, k is not small. For example, near the separatrix (see Fig. 1) it 
approaches unity. We have solved the general case in terms of the eUiptic integrals. In 
particular, in the vicinity of the separatrix the solution may be expanded over the small 
parameter k' — \/l — k"^ . These results will be published elsewhere. 
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APPENDIX A: 
DERIVATION OF EQUATIONS (5.2) - (5.4) 



In this Appendix we shall compute the derivative d < sin^ 9 >/ dT^in ■ Angle 9 is the 
one between the angular momentum vector J and the plane determined by the body's 
minor- and middle-inertia axes, 1 and 2. In the case of an oblate body, this angle also 
remains virtually unchanged through a period of precession, and changes slowly through 
many cycles (Efroimsky and Lazarian 2000)^^. In the general case of a triaxial body, this 
angle is not preserved through a precession period, though after one period of wobble 
it returns (adiabatically) to its initial value. The word "adiabatically" is in order here 
because in the course of many cycles angle 9 slowly decreases and thus deviates from 
the exact periodicity. 

In the body frame ( ei, e2, ) associated with the principal axes (1, 2, 3), the angular 
momentum components look: Ji = Qi Ii . Hence, according to 



|J| Sin^ = ^/2fi2 + _ ^J2^ ^ (J2 _ llS)ai . (Al) 

This shows that the angle 9 evolves with the same period as Vt\ = jS"^ sn'^[uj{t — to); k"^] . 
We remind that the periods of sn[uj{t — to); and cn[u}{t — to); k"^] are equal to 



4K{k^) = 4 (1 - k' sin2^)-i/2 = 27r (^1 + + 0{k' 



(A2) 



which is twice the period of dn[u!{t — to); k'^] . (See formula (16.1.1) in (Abramovitz and 
Stegun 1964)2^.) The A;^— approximation in the right-hand side of (A2) follows from the 
expansion (A7) below. The periods of the squares of sn and cn are equal to 



r 



2 2 K{e) =2 /"^' (1 - k^ sin2^)"i/2 ^^3) 







which is easy to understand from Fig. 16.1 in (Abramovitz and Stegun 1964) 
Squaring and averaging of (Al) yields: 

J2 < sin^ 9>= llR + {II - II S) <nl> = 

= i^hTurn - J') + h h ^ < sn\ > . (A4) 

— li ^ ' — li 

where u = uj {t — to) . The trajectories described by the angular velocity vector on 
the surface of the constant- J ellipsoid (Fig.l) are cyclic. The averaging may be performed 
over one such cycle, i.e., over the period 4K given by (A2). We shall though average 
over quarter of the total period, i.e., over K. This will be sufficient since in the above 
expressions (Al) and (A4) we have only squared components of f2 . Thus, 

< sn\ > = ^ (A5) 
K 

where the function Sn is, by definition (formula (16.25.1) in Abramovitz and Stegun 
1964)2^, squared sn integrated from zero to K. According to formula (16.26.1) in 
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(Abramovitz and Stegun 1964)27, Sn{K) = {- E{K) + K)/P where E{K{k^)) = 
/o"/' (1 - sin^ ey/'^de. Hence, 



2 1 EiK)\ 



(A6) 



Series expansions (formulae (17.3.11) and (17.3.12) in (Abramovitz and Stegun 1964)^'' 

(A7) 



K(P) 



TT 

2 



and 



TV 



1 • 3\2 



1 ■ 3 ■ 5 



2/ V2 ■ 47 3 V2 ■ 4 ■ 67 5 



entail: 



so that 



2 16 32 ^ ' 



y < sin' e> = ifR + {1^ - If s) < fi^ > 



(A8) 



(A9) 



^ (2/3 T.„ - J V ^2 ^3 1^ Hi + ^ + ^ + Oik^)) . (AlO) 



It follows from (AlO) and ( p.l4| ) that in the vicinity of pole C 



J < Sm y>(„earC7)— o " 77 TTTT U-'S^fcin — J I + 



+ T77 



16 (Is - /i)(/3 - hf y -2hT, 



+ 737 



(2 /3 Tfcj„ — J 



2\3 



kin 



32 (/3 - /i)(/3 - hf (J2 - 2/iTfc,„)2 



wherefrom 



d < sin 9 > 



+ Oik"") 



(All) 



dT, 



kin 



{near C) 



(J3 - Ji) (/3 - h) 



+ 



+ 



1 2 J3 Tkin — f I2 — II 



2 r2 



+ O(fc^) 



4 J2 - 2/irfe,„ V/3 - /27 /3 - h 
Together, the former and the latter yield: 

. fd <smH>\ h{hh + hh-2hh) 



(A12) 



dThi 



kin 



+ 



1 (/2 - hf 

2 (/3 - hfh - h) hh - hh - 2{h-h)h <sin'e> 



' (near C) 

/| < sin^ 9 > 



ih - h) ih - h) 

+ Oik") . 



+ 



(A13) 
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Now we shall derive similar formulae for the vicinity of pole A. Plugging (|2.16| ) into (AlO) 
we get: 

J2 < Sin^^ >inearA) = , , hT^in - j") + ^ , , (j' " 2 I.Tk 

I3 — II ^ ' Z I3 — Ii ^ 

+ 16 (2 J3r,,„ - J2) (/3 - - I,) + 32 (/3 - h)ih - hY (2 /3T,,„ - ^ ^ ^ 

This enables us to write down the derivative we need: 

2 f d < Sin^g > \ ^ /1/3 _ 1 hh{h - h) - 2hTuin . 4^ 

V dTk^n l^r^earA) ^ - h 4 (/3 - h) {h " h) 2 hTkin - ^ 

In the limit of oblate ( /3 > /2 = /i ) or prolate [I^^ = I2 > Ii) symmetry, expressions 
(A13) and (A15) simplify a lot: 



J 



7 -9/1 \ oblate T T / 1 • 9 n \ prolate ^ t 

2 I d < sm 6 >\ 2/3/1 t2 / < sm 6^ >\ h I3 



dTkin J (near C) \ ^Tkin J {near A) 



In the general case of a triaxial rotator, it ensues from (A14) and (A15) that 

j2 < sin^ ^ ^\ 

'''' ' +Oik^) . (A17) 



h - h 2 h - h h - h 2 < sin^ ^ > - 1 



APPENDIX B: 
SEVERAL WORDS ON THE QUALITY FACTOR 

In some situations it is possible to calculate the quality factor Q{uj) exactly. These 
are situations when one particular mechanism of attenuation dominates the others. For 
example, Q may be derived analytically for sound dissipation in a viscous liquid. It is said 
in (Landau and Lifshitz 1970)^^ that the calculation of the quality factor in a solid body 
would basically follow the same steps as in the case of a viscous liquid, in faith whereof 
the authors even present some laborious thermodynamical calculations. Unfortunately, 
in many cases this is not true, and the viscousity of solids contributes almost nothing 
to the attenuation (unless the body is warmed up to a plastic state). A much larger 
contribution to the attenuation is brought, in many materials, by phonon scattering on 
defects, and by a whole variety of related quantum effects (Nowick and Berry 1972)^°. 
In rocks, the attenuation is determined predominantly by displacement of defects. The 
numerous phenomena participating in the attenuation are subtle and bear a complicated 
dependence upon the temperature and frequency. Also mind a dramatic dependence of 
Q upon the humidity, as well as upon the presence of some other saturants. In many 
minerals, including for example silicate rocks, several monolayers of water may decrease 
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Q by a factor of about 55 (Tittman, Ahlberg, and Curnow 1976)^''. It is for this reason 
that the moonquakes cause an echo that keeps propagating and reflecting for long, almost 
without any attenuation. The echo would be dumped much faster, should the lunar 
lithosphere contain even a tiny fraction of water. Presumably, the moisture affects the 
inter-grain interactions in minerals. Another factor influencing Q is the confining pressure, 
but the pressure dependence is very weak within a broad (several orders) interval of 
pressures, and may be neglected. 

Returning to the frequency dependence, we would say that, fortunately, the overall 
frequency-dependence of Q is normally very smooth and slow, like for example, in the case 
of geological materials. Here follows the empirical temperature- and frequency- dependence 
of the quality factor, well supported by a vast experimental evidence (Karato 1998)^®: 

g ~ [cu exp(A7/?T)]" , (Bl) 

where A* may vary from 150 - 200 kJ/mol (for dunite and polycristalline forsterite) up to 
450 kJ/mol (for olivine). This interconnection between the frequency- and temperature- 
dependences tells us that whenever we lack a pronounced frequency-dependence, the 
temperature-dependence is absent too. At room temperature and pressure, at low fre- 
quencies (10~^ — 1 Hz) the shear Q— factor is frequency-independent for granites, and 
almost frequency-independent (except some specific peak of attenuation, that makes Q 
increase twice) for basalts (Brennan 1981)^^. It means that within this range of frequen- 
cies a is close to zero, and Q may be assumed also temperature- independent. For higher 
frequencies {lOHz — 1 MHz), the power a is remarkably frequency-insensitive (and equals 
approximately 0.25 for most silicate rocks). For a recent discussion and references on the 
Q-factor of asteroids see (Efroimsky and Lazarian 2000)^^. As for the Q-factor of the 
cometary nuclei, its value is unknown. Presumably, it should be close to the values of the 
Q-factor that are typical for snow and firn: between 0.5 and 5. 



APPENDIX C: 

AVERAGED OVER THE PRECESSION PERIOD SQUARES OF THE 
COMPONENTS OF THE STRESS TENSOR, IN THE VICINITY OF POLE C 



Formulae ( |8.7] ) - ( |8.17| ) trivially yield: 



<al>= ^ {x' - (l-Q)' Si 



^ (2I,T,,^ - J^)' l ^' - - ^'A ' (x' - a'Y h + 0{k')} 

32 V 3 ; I ^^^^ _ (j^ _ j^) J V ; I ^ 



PI 
4 



32 



2hTi 



kin 



h {h 



h (/l 



h h ih - h) ih - h) 



{y' - b'f {i + 0{k')} 
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<al>= I - c')' - Sf = (C5) 



2 P 

< at. > - 



' {(/37 + cxuk^) (y' - b') + (/37 - acvk') (x' - a')}' = (C7) 



xy ' ^ 

,2 /'o r T^. t2n2 



^ /./2/|(/3-/0(/3-/2) U^3 + /l-/2) -a)+ (a 

+ (/3 + /2 - /i) (l/^ - h')]' {l + 0(A;^)} , 



+ a7) (^^ - c^) + (-/3a; + a^) (x^ - a^)} S3 = (C9) 

(2/:-iI/,,:„, - J-) (J- - 2/i71.,„) r / 2 2\/r , 7- , / 2 2\/r , 7- r^\^ 

8" hllhih - h? -c){h + h-h) + [x -a){h + h-h)\ - 

-8¥/,/|/3(/3-/0M/3-/2)^^' -c)(/3 + /2-/0 + 

+ (x^ - a^) ih + I2- h)}'] {1 + 0{k')} , (CIO) 



{(a/3 + a;7) {z' - c') + {a (3 - u^) (y^ - 6^)}' 5, ^ (Cll) 



>2 (2/3Tfc,„ - J2) (J2 - 2/iT,, 



[8 /2 /a (/a - /i) (/a - /2) + - c ) + (/,- /a + /2)(y -6)} - 

Ilhh{h-h) {h-h? U^^ + ^3-/2)(. -c) + 

+ (/i - /3 + - {1 + 0{k')] , (C12) 

<(Tra)2>= [{x^ - a") {Q - I) + (y' - 6^) {Q + S) + - c^) (1-^)}' Si = (C13) 

.2 /o r_ -r2N2 
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[h h {h - h) - h h {h - h)] iy' -b') + [h h ih - h) 



hh {h - h)](z^-c')Y {l + 0{k')] . (CM) 



In the above formulae, factors Si , S2 , S3 , and S4 stand for averaged powers of the 
elhptic functions: 

Si = < (sn^(u, k^) - < sn^{u, /c^) > = 
= < sn\u, k^) > - < sn\u, >2 = 1 + 0{k^) , (C15) 



2 

> ) > 



S2 = < (^sn(u, /c^) cn(u, A;^) — < sn(u, fc^) cn(u, fc^) >^ 

= < sn\u, cn^[u, k^) > - < sn{u, k^) cn{u, k^) >2 = 1 + 0{k^) , (C16) 

8 

53 = < cn(u, k"^) dn(u, k^) — < cn{u, k"^) dn{u, k"^) >^ > = 

= < cn^(M, k'^) dn^{u, k^) > - < cn{u, k^) dn{u, >^ = (C17) 

54 = < (^sn{u, k^) dn{u, k^) — < sn{u, k^) dn{u, k^) >^ > = 

= < sn^{u, k^) dn'^iu, k^) > - < sn{u, k^) dn{u, >^ = (C18) 

where the averaging imphes: 

<...>=— / . . . du , 

r 7o 

r being the period expressed by (A2). The approximations were obtained by the trick 
( ^■181) - (|8.20|) explained in Section VIII. (Expansions of Sj over fc" cannot be obtained 
by plugging (|8.1| ) - (|8.3| ) into (C15) - (C18) because this would produce secular terms.) 

Plugging of (C2), (C4), (C6), (C8), (CIO), (C12) and (CM) into Q will lead us to 
the following expression for dissipation per unit volume: 
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d <W> d < 1^(2..) > ^ < ^(c.) > 

- + - 



dV 



dV 



dV 



(C19) 



where the first term stands for the dissipation associated with oscillations at the second 
mode: 

d < > 



dV 



- (or T - t2^' ^ 



1 + z/- 



- a^j (/i/2(/l - J2) - hh{h - h)) + [V' - (/2/3(/2 - h) - hhih - /l)) + 

+ [z^ - c^) (/3/i(/3 - h) - hhih - l2))f + {x^ - [hhih - h) - hhih - h)f + 

+ [y^ - b^y [hhih - h) - hhih - h)f + {z^ - c^y [hhih - h) - hhih - h)? + 

+ h h ih - h) ih - h) [ih + h - h) [x^ - a^) + (/3 + h - h) [y^ - 6')] '} + Oik^) (C20) 
while the second term stands for the dissipation at the frequency of precession: 

d < ly(^) > 



dV 



4/i 

~ 2 Ji Tkin) 

+ 



< > + < <^lz > 



^ J 2ij, 8 hhhih - h) ^ ^ ^ 



[hih - h] 
1 

h ih - h) 



ih - h)iz' -x' + a' - c') + h iz" + x^ - -c^ 



+ 



ih-h)iz' -y' + h' -c') + hiz' + y' 



1 p2 



r2\2 h - h 



1 



+ 



h - h [hih - h 
1 5 



2// 8 hhhih - h) 
ih - h)iz^ - + - c^) + h iz^ + - a 



(2/3T; 



kin 



2 C^l 



+ 



hih 



ih-h)iz'~y' + b''-c') + hiz' + y' 



+ Oik^) . (C21) 



As expected, we have obtained the result, (C20 - C22), in the spectral form ( |6.8|) . After 
integration (like (p.9|)) this result must be plugged into ( |5.4| ), which will lead to ( |8.24| ). It 
follows from (C20) and (C21) that the geometrical factors emerging in ( ^.241) will read: 
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1 1 f"^ ( 1 
Ht = ; / dx dy dz { — ; 

8 Ji /2 h {h - h) J -a J-h J-c \ h {h - h) 



ih - h){z^ - + 0^ - 0^) + 



I r ^ 2 I 2 2 2 

+ hyz + X — a — c 



+ 



+ 



1 



h ih - h 



Ho 



1 r ( 1 

dx dy dz { 

J-b ^ J-c I 1 + 



(C22) 



X 



' 64 /f/|/|(/3 - h)Hh - h? 
- a^) {hhih - h) - hhih - h)) + {y^ - h^) {hh{h - h) - hh{h - h)) + 



+ (z^ - c^) (/3/i(/3 - h) - hHh - h))] + {x^ - a^) [hhih - h) - hhih - h)? + 
+ (y^ - b'Y [hhih - h) - hhih - h)? + {z^ - [hhih - h) - hhih - h)? + 



+ hh ih - h) ih - h) [ih + h - h) (x' - a") + (Is + h - h) {y^ - 
and 

3 



1 2 



(C23) 



1 /1/3 + 12/3 - 2/1/2 h - h r, f' r , J 

64 ih-h)Hh-hy hhh J-a'^'^J-byj-cXhih-h) 



c^)ih + h - h) + ix' - a^)ih + h- h 



+ 



h ih - h 



i^' 



c')ih + /a - h) + iy' - b')ih + h- h) 



For a homogeneous prism of dimensions 2a x 2b x 2c : 



a 6 



6^ a' 
+ 



(62 + c2) (a4 - c4) (a^ + b^) \b^ - - 



(C24) 



(C25) 



100 a-^b^c - aH^c^ - aH'-^c^ + 0.21 aHc^ + 0.19 aH^c^ + 0.21 a b^c^ 
-"2 - — r , o ,os2 . . 772 (C26) 



and 



Ho 



237 



T 5/2 

a b c [a 



(2.67aH4 - a^c^ - l.Q7b^ c'^) ia'^ + 6^ - 2c 



(C27) 



• ' (a2 + 62) (^2 _ c2) (62 - c2) (a4 - c^y (6^ - 0^)2 ' 

Calculating Hi and //2 we assumed that the quantity 1/ (1 + z/^^) emerging in ( |7.2| ) and 



S.22|) is 1/5 (as it normally is for solid materials). We also used the standard formulae 
for the moments of inertia: 

A = y (6' + c' 
m being the mass of the homogeneous body. 



etc... 



(C28) 
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APPENDIX D: 

AVERAGED OVER THE PRECESSION PERIOD SQUARES OF THE 
COMPONENTS OF THE STRESS TENSOR, IN THE VICINITY OF POLE A 



By squaring each of the expressions 
easily arrive to the following formulae: 



( p.l4| ), and averaging the result, one will 



< o^^ > = — 



kin I 



32 II II (Ji - I,Y (I, - I,) 



{I, ih - h) - h (/i - h)f {x^ - a^Y {l + 0{k')] , 



32 72/32 (/, - i^f - i^f ^ ' ^ ' ' ^ ' ) \ ^ ^ 



z - C ^1 



(J2 - 2 Ji T, 



kin ) 



32 /2/2 (/^ _ i^y (/, 



-J^ {I2 ih - /2) - h ih - h)V - C') {1 + Oik')} 



^ ^2 ^ P 



(D7) 



(J^ — 2 /i Tfcj„) (2 /s Tfcj„ — 



j't P T UJ 7^ (x^ + y^-a^ - h-) + (D8) 

+ (/2 - h) {x' -y' + b' - a')Y {1 + Oik')} , 



/9 



(D9) 



(2^3T. - J^) (J^- 2 /,r..„) .^^^ _ _ ,2 + ^2 _ ^2) + 

[ O Ji i2 -'3 Us — -'ij ^ 
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1^3 — ^2) K'f — ^ J-kiri) (fj T \ ( 2 2,2 2\ , mi 



3 p2 _ J^) (J2 _ 2 Ji Tfc, 



32 



< {Traf >= P-P' {{x^- a') {Q - 1) + {y' - b') {Q + S) + {z' - c') (1 - S)}' 



(D13) 



(J2 - 2hTf 



kin ) 



32 /2 /| /| (Ji - (/3 



{[/l/2(/l-/2)-/l/3(/l-/3)] 



- ) + 



+ [hh ih - h) - hh ih - h)] [y' - h') + [/3/1 (/3 - h) - 

- hh {h - h)]{z^-c^)Y {1 + 0{k')]. (DM) 

where ^1^2,3,4 are given by (C15) - (C18). Just as in the preceding case of pole C, the 
above expressions (D2), (D4), (D6), (D8), (DIO), (D12), (D14) form pole A are to be 
plugged in ( [7l2|) . It will entail: 



d <W > d < > d < > 



dV 



dV 



dV 



(D15) 



where 



dv = IJll^ +2 



1 (J^ — 2 h Tkin) (2 h Tkin — J^) 

271 "s" 



1 



hhh{h - h] 



h{h - h 



ih - h) {x^ - - + c^) + 
+ h [x'^ - + z'^ - < 



+ 
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+ 



and 



+ h (x^ - + - c^)] + Oik"") (D16) 



1 

4/i 



< [Traf > + < > + < > + < f^L > +2 < crl > 



yz 



1 p2 ( j2 _ 2 Ji 



4/i 32 /2/|/|(/3-/i)2(/2-/i 



[/i/3(/i - Is) - hUh - I2)? (x^-a^y + 



+ [hh [h - h) - hh {h - h)f iy" + 



+ [13/2 (/3 - h) - hh ih - h)? (z' -C'V + 



+ 2hh{h-h){h-h) [{h + h + h)iz-c') + ih-h + h){y'-b') 
- Y^—T {[hhih- h) - hhih- h)]{x'- a') + [hh ih - h) - 
h h ih - h)] {y' - b') + [h h ih - h) - h h ih - h)] {z' - c') Y] + Oik') . (D17) 



Integration of the two above expressions over the volume of the body gives expressions for 
the two components of the time- dependent elastic energy deposited in the body: W^"^^ 
and W^"^^^ , plugging whereof into (|6.6| ) will yield ( |9.15|) . Expressions for the geometrical 
factors emerging in ( p.l5| ) are: 



1 1 f"^ f'' i 1 

16 hhhih - h) J-a'^'' J-b'^^ J-c'^^ [hih - h] 



Si 



ih - h) [x' - a' - + c'] + 



+ h [x^ - + - c^] ' + 



(D18) 



+ 



h ih - h 



h {x' -a' + y' - b') + ih - h) {x' -a'-y' + b') 



(D19) 
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+ [/2/i(/2 - h) - hhih - h)f {y^ - b'Y + 

+ [13/2 (/3 - h) - hh {h - h)? (z^ - C')' + 

+ 2/2/3 (/3 - h) ih - h) [ih + I2 + h) {z" - c') + (/i - /3 + h) {y' - b')]' - 

- {[^1^2 (/i - h) - hh {h - h)] {x' - a') + 

+ 1/2/3 ih - h) - hh ih - h)] {y' - b') + [/3/1 (/3 - h) - hh ih - h 



The factor 5*0 is equal to 



A/|/,(W.n/.-/.) I> l> I> {(^» - i^' -"'-^'^ + 

+ h{x' - a" + - c^)Y (D20) 



ra fb 



128 



and becomes negligibly small in the case of the body approaching the prolate symmetry 

(/3 - /2 ^ 0). 

For a homogeneous prism of sizes 2a x 2b x 2c , the factors much simplify: 



^1 



173 



(I b ( 



(a2 + 62) (^,2 + (a^ - c^) 



+ 1.7 « ^// + If + 1.7 (fc^ + c' 



a4 - 64 



+ 



a4 - c4 



, (D21) 



ah c 



+ 2.8 a^^c^ - 4.80^^0^ + 

4 l4 



- 4.8 aH* c« + 4.8 6' c« (D22) 



and 



32^ (gS + 1.67 g^c^ + c') (62 - c2) 
° ~ (g4 - 64) (g4 - c4)2 (62 + c2) 



(D23) 
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